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Abstract 

We show how generalized quark distributions in the nucleon describe the density of polarized 
quarks in the impact parameter plane, both for longitudinal and transverse polarization of 
the quark and the nucleon. This density representation entails positivity bounds including 
chiral-odd distributions, which tighten the known bounds in the chiral-even sector. Using 
the quark equations of motion, we derive relations between the moments of chiral-odd 
generalized parton distributions of twist two and twist three. We exhibit the analogy 
between polarized quark distributions in impact parameter space and transverse momentum 
dependent distribution functions. 



1 Introduction 

The distribution of transverse quark spin in the proton remains one of the most intriguing and least 
known aspects of nucleon structure. It has been the subject of numerous theoretical studies, and 
there is a vigorous experimental program aiming to measure the transversity distribution h\(x) in 
present or planned experiments. Recent overviews and references can for instance be found in pQ. 
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A wealth of information on the nucleon structure is encoded in generalized parton distributions 
(GPDs), see e.g. the reviews [21 El ID- They admit a particularly intuitive physical interpretation at 
zero skewness £, where after a Fourier transform they describe how partons with given longitudinal 
momentum are spatially distributed in the transverse plane pj- A remarkable spin effect in this 
representation is that transverse nucleon polarization induces a sideways shift in the quark density, 
whose size is related to the anomalous magnetic moment of the nucleon and thus quite substantial [5] . 

A relatively small number of studies have so far been devoted to generalized transversity distribu- 
tions, which were introduced in [71 012]. Since the operator measuring transversity is chiral-odd, it is 
notoriously difficult to find processes where transversity distributions can be accessed experimentally. 
For generalized transversity distributions it is indeed not clear if this can be achieved in practice, and 
at present there is only one type of process known where this may be possible in principle |1U) . There 
is however the prospect of gaining information from lattice QCD, which provides a tool to calculate 
the Mellin moments of generalized parton distributions. Several studies have been performed for 
chiral-even distributions |11|I12|. and first results for chiral-odd ones have been presented in The 
purpose of this paper is to take a closer look at the physical interpretation and properties of these 
quantities. 

In Sect. [21 we extend the analysis of [S] to generalized transversity distributions and investigate 
the distribution of transverse quark polarization in the impact parameter plane. The result closely 
resembles the expressions for the distribution of polarized quarks as a function of their transverse 
momentum. In Sect. El we derive positivity bounds which involve chiral-even and chiral-odd distri- 
butions, both in impact parameter and in momentum representation. We also give bounds that are 
valid for Mellin moments. Sect. [U is devoted to relations between distributions of twist two and three 
resulting from the QCD equations of motion. We summarize our findings in Sect. El 



2 Polarized parton distributions in the transverse plane 

To begin with, let us recall the definitions for generalized quark distributions in the proton. Following 
the conventions of |14[ El E] the distributions of twist two read 
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Corresponding to the quark-antiquark operator in their definition, the distributions parameterizing 
F and F are referred to as chiral-even, and those parameterizing F^ as chiral-odd. The latter are 
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also called quark helicity flip or generalized transversity distributions. We use light-cone coordinates 
= (v° ±f 3 )/\/2 for any four- vector v and write its transverse part as v = (v 1 , v 2 ). Scalar products 
of boldface vectors are defined such that v 2 > 0, and Roman indices i, j, k are understood to be 
restricted to the two transverse directions. Our sign convention for the totally antisymmetric tensor is 
^0123 = 1- We use kinematical variables P = \{p + p'), A = p'—p, t = A 2 , £ = (p — p') + /{p+p') + and 
denote the proton mass by m. For better legibility we have not explicitly labeled the polarization of 
the proton states {p'\ and \p) and have omitted the momentum and polarization labels of the proton 
spinors u and u. The definitions Q hold in the light-cone gauge A + = 0, otherwise a Wilson line 
appears between the quark field and its conjugate. 

We will find that in all expressions of this paper the distribution Et appears in the combination 
Et + 2Ht, so that one may regard Et + 2Ht as a more fundamental quantity than Et- Using the 
Gordon identity one can make this combination appear already in the decomposition of the matrix 
element F T (x,£,t), rewriting 



Ht ua +J 75 u + Ht u 2 - u + Et u o _ - — u (2) 
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In this and the next section we restrict ourselves to skewness £ = 0, where generalized parton 
distributions have a probability interpretation when transformed to impact parameter space jSj. To 
make this explicit we form wave packets 

|p+ ' b)= /(2^ 6 ^ P|p) (3) 

from the states states \p) with definite four-momentum, where it is understood that the integration 
over p is done at fixed p + with p~ = (m 2 +p 2 )/{2p + ). The state \p + , b) has definite plus-momentum 
p + and definite impact parameter b, i.e., it is localized at position b in the x-y plane. Further analysis 
shows that b is the "center of momentum" of the partons in the proton |15j . given in terms of their 
plus-momenta and transverse positions as b = J2ipfb { / J2iPt • A two-dimensional Fourier transform 
gives 

F(x,b) = [ ^- e - ibA F(x,0,-A 2 ) 



(2nf 



-i 



/ ^ elXP+Z ~ (p + ,0\q(z2h + q(z 1 )\p + ,0) (4) 



with zf = = 0, z\ = Z2 = b, and = — z^ = \z~ . Here we have used translation invariance to 
shift the quark-antiquark operator to transverse position b and the impact parameter of the proton 
to the origin. The normalization factor J\f = (2ir)~ 2 j d 2 p is singular like a delta- function, which can 
be avoided if instead of © one takes wave packets smeared out in impact parameter space I n 
analogy to (@J) we define matrix elements F(x,b) and F T (x,b). The impact parameter distribution 

H(x, b 2 ) = f e- ibA H(x, 0, - A 2 ) (5) 



(2ny 

and its analogs E(x, b 2 ), H(x, b 2 ), Ht(x, b 2 ), etc. depend on b only via its square thanks to rotation 
invariance. We see in (jlj) that the Fourier transformation has made the matrix element diagonal in the 
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plus-momentum and the impact parameter of the proton states. If we also take the same polarization 
for these states, the matrix element becomes an expectation value and thus acquires a probability 
interpretation akin to the usual parton densities. 

The wave packets (J5J involve proton momenta which are not along the z-axis, and the spin 
states for this case have to be chosen with some care. It is useful to take states of definite light-cone 
helicity ^7]. A proton state of momentum p with positive (negative) light-cone helicity is transformed 
to a proton state at rest with spin along the positive (negative) z-axis by a Lorentz transformation 
C(p), which is the combination of a transverse and a longitudinal boost (see Sect. 3.5.1 of [3] for a brief 
summary). The light-cone helicity of a state is invariant under boosts along the z-axis, and for large 
p + light-cone helicity coincides with the usual helicity up to effects of order m/p + . The superposition 
( |+) + e lip \— ) )/\/2 of states with positive and negative light-cone helicity is called a state of definite 
transversity, which can be seen as the light-cone analog of definite transverse polarization. According 
to what we just discussed, C(p) indeed transforms this state to a state at rest whose spin vector is 
given by S = (cos (p, sin ip) and S z = 0. A state with both longitudinal and transverse polarization 
can be written as 

|A, S) = cos(±tf) |+) + sin(±i?) e lLp \-) (6) 

and is transformed by C{p) to a state at rest with spin vector given by S = (sini? cos(/?,sin?9 sin</?) 
and S z = cost?. We will therefore use S and A = S z to characterize these states. Combining them to 
wave packets (0) we finally obtain states suitable for interpreting the matrix elements F(x, b), F(x, b) 
and Fj,(x, b). For ease of language we will call S and A the transverse and longitudinal polarization 
of the proton. In the following it will be important that they respectively transform like a usual spin 
vector and usual helicity under rotations in the x-y plane and under parity or time reversal. 

For quarks and antiquarks we consider light-cone helicity states as well. Note that the quark 
operators in (@J are at definite transverse position and thus correspond to integrals over the quark 
or antiquark transverse momentum. Quarks with light-cone helicity A = ±1 are projected out by 
the operator ^7 + [l + A75] q. Evaluating the proton spinor products in £Q) for the states © and 
Fourier transforming the result, we obtain the density of quarks with light-cone helicity A, light-cone 
momentum fraction x and transverse distance b from the proton center as 

1 



1 



F(x,b) + XF(x,b) 



H(x, b 2 ) - S' l e ij V- -J^E{x, b 2 ) + XKH{x, b 2 ) 
m ob A 



(7) 



for x > 0, where repeated Roman indices are to be summed over. For x < the density of antiquarks 
with light-cone helicity A, light-cone momentum fraction — x and transverse position b is given by 
\[— F(x,b) + \F(x,b)]. It readily follows from the transformation properties of qj + q and g7 + 7s q 
under charge conjugation that in going from quark to antiquark densities one has to change the sign 
of F but not of F. The result Q is well-known and for instance discussed in The term with H 
reflects the difference in density of quarks with helicity equal or opposite to the proton helicity. More 
remarkably, the term with E describes a sideways shift in the unpolarized quark density when the 
proton is transversely polarized. 

We now discuss transverse quark and antiquark polarization, which in analogy to our above dis- 
cussion we define as the superposition ( | +) + e lx \ — ) ) /\/2 of positive and negative light-cone helicities, 
with a transverse spin vector s = (cosx>sinx)- Quarks with transverse polarization s are projected 
out by the operator ^g7 + [l + (57)75] q = \q [7 + — sHa^j^] q, and their density is 



F + s l F\ 



H - SVV- E' - sW- ( E' T + 2H' T ) 



+ s 



sHht 



Am 2 



A b H T ) +s i (2b i b j 



(8) 
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Figure 1: Density plots in the impact parameter plane for the functions 6exp[— b 2 /b^] sin eft (left) and 
6 2 exp[— 6 2 /6q] cos(20) (right), with bo = 0.5 fm. These functions illustrate the terms in the quark 
density (|SJ) which break rotational symmetry, as explained after Dark areas represent high 

densities. 

for x > 0. The density of antiquarks with transverse polarization s and light-cone momentum fraction 
—x is given by — ^[F(x, b) + s l F^(x, b)]. Here and in the following it is understood that when nothing 
else is indicated, the matrix elements F, F, Fj, and the distributions H, E, H, Ht etc. are functions 
of x and b as given in (j3J) and ©. We write b = \fh 2 so that b 2 = b 2 and d/db 2 = d/db 2 , and we use 
the shorthand 

for the derivatives and 

«4w'=4(4)' < 101 

for the two-dimensional Laplace operator acting on functions / that depend on b only via its square. 
In (jHJ) we have further introduced the two-dimensional antisymmetric tensor e 4 - 2 with e 12 = — e 21 = 1 
and e 11 = e 22 = 0. 

The term with E' T + 2H' T in (jHJ) describes a sideways shift in the distribution of transversely 
polarized quarks in an unpolarized proton, whereas the last two terms in that expression reflect 
a correlation in the quark density between the transverse polarizations of quark and proton. The 
structures which break rotational symmetry in the density (jHJ) are 

STftiV = bsm4>, 
sVV = bsm(<fi-x), 
s i (2b i V -b 2 5 ij )S j = 6 2 cos(x-20), (11) 

where we have parameterized b = b (cos (ft, sin (ft) and taken S = (1,0) for simplicity. For illustration we 
show in Fig.^density plots for the functions 6exp[— 6 2 /6q] sin</> and b 2 exp[— 6 2 /6q] cos(2</>), where the 
exponentials are taken to mimic the impact parameter dependence of the relevant parton distributions. 
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Integrating the densities (JJJ) and (JBJ) over all impact parameters one obtains generalized parton 
distributions in momentum space at t = 0, namely 



/ 



d 2 b 



d 2 b 



F(x,b) + XF(x,b) 
F(x,b) + s l F^(x,b) 



= H(x,0,0) + XAH(x,0,0) 
H(x,0,0) + s i S i H T (x,0,0) 



h(x) + XA gi (x) , (12) 
f 1 (x) + s i S i h 1 (x). (13) 



Here we have used that GPDs in the forward limit £ = and t = reduce to the usual parton 
densities, namely 



if (x, 0,0) = fi(x), 



H(x,0,0) = 9l (x), 



H T (x,0,0) = h^x) 



(14) 



for x > 0, where f\ denotes the unpolarized quark distribution, g\ the quark helicity distribution, and 
h± the quark transversity distribution (another common notation is f\ = q, g\ = Aq and hi = 5q). 
Corresponding relations involving antiquark distributions hold for x < 0. Weighting the impact 
parameter distributions with b 2 before integration, one obtains derivatives at t = 0, 



Jd 2 bb 2 (F + XF) = A^H + XAH 



d 2 bb 2 (F + s i F^ 



-( 



dt 



H + s l S* 



4m 2 



Hj 1 



(15) 
(16) 



J t=o 



where we use the subscript t = to indicate that the GPDs are taken in momentum space as in 
with £ = as always in this and the next section. The ratio of the integrals in (|15|) and (|12|) or in 
(|16|) and ()13|) thus gives the average squared impact parameter b 2 of quarks with given polarization 
and plus-momentum fraction. The average sideways shift in the impact parameter distribution due 
to the transverse polarization of either the quark or the proton is obtained from 



J d 2 b V (f 



+ s l F 



r 



— ( S^^E + sV J '(£ T + 2H T ) 
2m V 



(17) 



J t=o 



normalized to the integral in (|13|) . This shift is more involved than in the case of longitudinally 
polarized or unpolarized quarks, 



d 2 bb> [F + XF 



d 2 b V F 



2m 



(18) 



J t=o 



which has been discussed in some detail in [S] . Finally, the average distortion of the impact parameter 
density due to the last term in ® is characterized by 



J d 2 b (2b j b k - b 2 5 jk ) (F + s l F l T 



-At (s j S k + S j S k 

m 



s i 3 i S jk^ Ht 



(19) 



J t=o 



We note in Q and (jHJ) that there is no polarization effect in the impact parameter distributions 
for longitudinally polarized quarks in a transversely polarized proton and vice versa. This is because 
the only structures describing such effects which are allowed by parity conservation are XS l b t or 
As l b l . They are odd under time reversal and hence forbidden. This corresponds to the fact that 
the generalized distributions E and Et in (^Q) do not appear in the matrix elements at £ = 0, the 
former because it is multiplied with A + = — 2£P + in its definition and the latter because it is an odd 
function of £ (Hj • 
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It is instructive to compare our impact parameter densities with the densities for quarks of definite 
light-cone momentum fraction x and transverse momentum k, which play a prominent role in the 
description of spin asymmetries in a variety of hard processes. They can be defined from the correlation 
function 



dz d 2 z 
4vr (2vr) 2 



-ikz 



ip\qp{-\z)W[-\z,\z] q*{\z)\p) 



(20) 



where it is understood that the proton states have zero transverse momentum. The Wilson line W 
has recently been recognized as essential in the definition, since different physical processes require 
different paths leading from —\z to \z and can actually give different correlation functions (see e.g. 
|18j and references therein) . Physically, the gluons resummed in the Wilson lines describe interactions 
of spectator partons in the process where the correlation function appears, and the corresponding 
parton distributions describe the density of quarks or antiquarks in the presence of these interactions. 
This subtlety did not appear in the generalized parton distributions (J2) and their impact parameter 
analogs, because there the quark field and its conjugate are separated by a light-like distance and the 
relevant Wilson line just runs along the light-cone between —\z and \z. Projecting out densities for 
quarks of definite longitudinal or transverse polarization from Q2U|) . one obtains ^5] 
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(21) 



where we have used the notation of Boer, Mulders, Tangerman |2()1I21| for the distribution functions, 
which depend on x and k 2 . Integrating over k one recovers the distributions fi(x) = J d 2 k fi(x, k 2 ), 
gi(x) = j d 2 k gi(x,k 2 ) and h\(x) = J d 2 khi(x,k 2 ) we already encountered in (|14|). The tensor 
structures in ()21j) are analogs of those in Q and (jHJ), with k taking the role of b. The corresponding 
analogy between transverse momentum dependent and impact parameter dependent distributions 
reads 



fx 



A b H T /(4m 2 



/it 
hi 



E', 

■ i^E'rp ~\- 2Hrp) 



9i « 
n 1T 



* 2HZ 



(22) 



The impact parameter distributions which would correspond to g lT and h^ L vanish because of time 
invariance, as discussed above. Notice that the momentum k changes sign under time reversal, 
whereas the position vector b does not. The structures \S l k l and Ks l k l describing polarization 
effects for longitudinally polarized quarks in a transversely polarized proton and vice versa are hence 
time reversal invariant. On the other hand, both S l e l3 k 3 and s % e ll k? are odd under time reversal. The 
corresponding distributions f^t an( ^ ( w hich respectively are the Sivers and Boer- Mulders functions) 
are however not constrained to be zero by time reversal invariance. This is because the Wilson 
line in the correlation function for relevant processes like semi-inclusive deep inelastic scattering or 
Drell-Yan pair production have paths that are not invariant under time reversal, contrary to the 
paths appearing in the impact parameter distributions discussed above. Time reversal thus connects 
transverse momentum dependent distributions with different Wilson lines, but does not constrain 
them to be zero |22| . We finally note that, beyond the formal correspondence of the functions E(x, b) 
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in iJBJ) and f^p(x,k) in (|21jl . a deep dynamical connection between them has recently been proposed 

in HH1I231 

We have seen in (|12j) to (|19|) how the impact parameter distributions can be reduced to distribu- 
tions only depending on the momentum fraction x by taking appropriate integrals over b. Conversely, 
one obtains distributions only depending on b by integrating over x. Taking Mellin moments in x, 
we obtain expectation values of the well-known local twist-two operators in proton states localized at 
zero impact parameter, 



{p + ) n J 1 dxx n ~ l F(x,b) = ^AT 1 (p+,o|g 7 + (iB+) n - 1 g|p+,0^ 

{p + ) n J l ^dxx n - l F{x,b) = i AT 1 | ^75 {iD + ) n -' ■ 

(p + ) n J^dxx n - l F^x,b) = -^AT 1 (p+,0\qa +:i l5 (iD + ) n - 1 q\p + ,0), (23) 

where = \{D^ — D^) = \{d^ — d^) — igA^ , and all field operators are to be taken at position z 
with z + = z~ = and z = b. To obtain matrix elements of local operators, one has to integrate over 
both positive and negative x and hence not only combines the information from different momentum 
fractions but also from quarks and antiquarks. According to the charge conjugation properties we 
discussed after and (jHJ), moments with odd n in (|23j) correspond to the sum of quark and antiquark 
densities for F and to their difference for F and FiL, whereas for moments with even n the situation 
is reversed. The lowest x moments of F and ^ [F + s l F^] hence describe the transverse distribution of 
unpolarized and transversely polarized quarks minus antiquarks in the proton, respectively. Higher x 
moments give the transverse distributions of quarks plus or minus antiquarks weighted with a power 
of their plus-momentum fraction. In contrast, the x moments of ^[F + \F] describe the transverse 
distribution of quarks plus or minus antiquarks with chirality A (i.e. quarks with helicity A and 
antiquarks with helicity —A). A two-dimensional Fourier transform turns the expectation values (j23|) 
into matrix elements for proton states of definite momenta, which are parameterized by form factors 
depending on the squared momentum transfer t. These form factors can be evaluated in lattice QCD 
since the corresponding operators are local and thus allow continuation into Euclidean space. It is 
amusing that after a Fourier transform they become quantities (|23[) whose physical interpretation is 
naturally given in a light-cone framework. 

Note that our interpretation of form factors differs from the well-known interpretation due to 
Sachs j2E]) where their three-dimensional Fourier transforms yield densities in a static proton state. 
That framework has been extended to generalized parton distributions as functions of x, £ and t in |26| . 
It is limited by ambiguities due to the impossibility to localize a particle more accurately than within 
its Compton wavelength. Localization in only two dimensions is not affected by this limitation, and 
the wave packets © are indeed eigenstates of a two-dimensional position operator [T2|. The price to 
pay for this is the loss of manifest three-dimensional rotation invariance in the light-cone framework. 
In return, the mixed representation of position space in two dimensions and plus-momentum in the 
third allows one to boost to a frame where the proton moves fast, which is a natural frame for the 
physical interpretation of quark and antiquark degrees of freedom. 

So far we have interpreted ^[F + s l Fj] in (JHJ) as the density of quarks with a given transverse 
polarization s. The vector field -FY(x,b) gives the direction in which the transverse polarization of 
quarks is largest, and its size \FT{x,b)\ is the difference of densities with quarks polarized along or 
opposite to this direction. According to @ we can write F^ as the superposition of three terms, 
given by functions of b 2 times the vectors S z , e^lfl and {2b l lP — b 2 8 l ^)SK The field lines of the term 
with e^W are circles around the origin, and those of the term with (2b l b^ — b 2 5 % i)Si are circles going 
through the origin with a tangent along the proton polarization S 1 , as illustrated in Fig. [2j 
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Figure 2: The vector fields exp[-6 2 /^] (left) and (26*6* - b 2 5^)S^ exp[-6 2 /^] (right) with S l 
along the x-axis and bo = 0.5 fm. They illustrate the form of two terms in the decomposition of 
Fj,(x,b), which describes the transverse polarization of quarks in the impact parameter plane. The 
third term in the decomposition is a field parallel to S l . 



The divergence and the curl of the field F^(x, b) respectively are 




:R 

db* T 



2S l b l H' T , 



(24) 



They can be rewritten as matrix elements of operators which are total derivatives. This is readily 
seen in Mellin space, where we have 



i 



--M- L {p + ,0 d fM [qa + ^ 5 (zD + 
AT 1 (p+, SL Iqa^iiD+^q 



\n-l. 



q 



p+,0 



p + ,o , 



(25) 



with all field operators evaluated at z + = z~ = and z = b. To obtain (|25[) we have used the 
representation (|2U)) and the fact that the first term in d + [qa^ . . . q] + d{ [qa +l . . . q] = [qa +>1 . . . q] 
vanishes when taking a matrix element between states of equal plus-momentum. The operators in 
(|25|) can be rewritten using the equations of motion as we will discuss in Sect.0J 

The Mellin moments of the impact parameter distributions H(x,b 2 ), E(x,b 2 ) etc. are Fourier 
transforms of form factors in momentum space, which are denoted by 

A n0 (t) = J dxx n ~ 1 H{x,Q,t), B n0 (t) = J dxx n ~ 1 E(x,G,t), 

A Tn0 (t) = J dxx^HrfaQJ), B Tn0 (t) = J dajz n-1 £>r(x,0,t), 



A Tn o(t) 



dxx n - l H T (x,0,t) 



(26) 
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in a standard notation (see Sect.EJ). These form factors can be calculated in lattice QCD ^DEIEI]; 
where it has become customary to fit them to a dipole form. For reasons that will be clear shortly, 
let us consider the more general power-law ansatz 

A{t) = n f/% ' (27) 

(1 - t/m 2 A Y 

where the power p and the mass m A are free parameters for a given form factor A(t). Note that in 
the limit p — » oo at fixed m\/p this ansatz gives an exponential in t. The Fourier transformation of 
(|2Tj) to impact parameter space leads to the modified Bessel function, 

A(b 2 ) = C (mAbf^K^imAb) , C = A(0) , (28) 



and the derivatives defined in © and Q1U|) are 

A'(b 2 ) = -\Cm 2 A {m A by- 2 K p „ 2 (m A b), A" \b 2 ) = \Cm\{m A bf^K p ^(m A b) , 

A b A(b 2 ) = -Cm 2 A (m A b) p - 2 [2K p ^ 2 {m A b) - m A b K p ^ 3 (m A b)j . (29) 

A parameterization of the type (|27|) is in the first instance only valid in the t range where the form 
factor has been fitted. In particular, lattice computations have an upper limit |t| m ax on the squared 
momentum transfer given by the lattice parameters. This corresponds to a limited resolution of order 
(l^lmax) -1 / 2 on the impact parameter [Ej- Furthermore, results obtained on a finite lattice cannot 
give direct information on the behavior of quark densities at impact parameters much larger than 
the lattice size. Nevertheless, one may want to require that a parameterization leads at least to a 
physically plausible behavior of the impact parameter density at small and large b. To analyze this 
behavior we need the relations 

2 2 p ~ 1 T(v) 
K (z) zZq log-, K p{z) forp>0, 

and K-p(z) = K p (z). At large b, each term in the Mellin moments of the densities Q and (JHJ) 
then falls off like (m A b) p ~ 3 / 2 e~ mAb . For the limit 6 — > it seems reasonable to require a regular 
behavior of the impact parameter density, which implies that no term should diverge and that WE 1 , 
IP(E' T + 2H' T ) and (2b 1 \P — b 2 5 l ^)H'^ should vanish at b = 0, since they have a nontrivial dependence 
on the azimuthal angle <j). This restricts p in the parameterization of moments to p > 1 for H, H and 
Ht, to p > 3/2 for E and Et, and to p > 2 for Ht- The terms with H, H, Ht and A^Ht in the 
moments of (J2J and (JHJ) then all take finite values at b = 0. In momentum space these restrictions are 
tantamount to requiring that the Mellin moments of F(x,0,t), F(x,0,t), Fj,(x,0,t) decrease faster 
than 1/t for t — ► — oo, as is readily seen when the proton spinor products in (QJ are evaluated Note 
in particular that a dipole ansatz with p = 2 for ^4x^0 (t) gives only a 1/t falloff in the nth moment of 
Fj,(x,0,t) and a corresponding logarithmic divergence at b = in the nth moment of Fj>(x,b). 

To illustrate how the transverse spin density in the proton may look like, we focus now on the 
first moment \ j\ dx [F(x, b) + s l F^(x, &)], which gives the difference of impact parameter densities 
for quarks and antiquarks (for ease of language we will simply speak of quarks in the following). As a 
numerical example we take a parameterization (|2*Tj) with p = 2 for Aio(t), Bio(t), Axioit) and p = 3 
for Bxioit), Ayio(i). We set the mass parameters m A to 1 GeV for Aio(t), Bio(t) and to 1.5 GeV for 
the three other form factors, and take 

^io(O) = 2, Bio(0) = 3, Ano(0) = 1, Bno(0) = 6, A T io(0) = -1 (31) 
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Figure 3: Left: Illustration of the first moment f_ 1 dx F(x,b) of the impact parameter density for 
unpolarized u-quarks in a proton with transverse spin vector S = (1,0). Right: The same for the 
first moment \ I-i dx [F(x, b) + s l F^(x, b)] of the distribution of u-quarks with transverse spin vector 
s = (1, 0) in an unpolarized proton. Dark areas represent the highest and light areas the lowest values 
of the density. Further explanation is given in the text. 

for their values at t = 0. This set of parameterizations is a rough approximation of preliminary 
results from lattice calculations [27j for the first moments of generalized u-quark distributions (where 
\t\ goes up to about 3.5 GeV 2 ~ (O.lfm)" 2 and lattice sizes are between 1.5 and 2.2 fm). We stress 
that it is meant to be indicative and not a precise representation of those results. We note that 
^4io(0) = 2 correctly gives the number of valence u-quarks in the proton, whereas -Bio(O) = 3 is too 
large compared with the value 1.67 one obtains from the measured magnetic moments of proton and 
neutron (recall that Bio(t) is the relevant quark flavor contribution to the electromagnetic Pauli form 
factor). 

In Fig.|3]we show the resulting first moment of the impact parameter density for unpolarized quarks 
in a transversely polarized proton and for transversely polarized quarks in an unpolarized proton. The 
dipole-type structures due to S l ^WF' and s l ^W(E' T + 2H' T ) are clearly visible, reflecting the large 
values in 1)3 lj) for B\q and Btio + 2Atio at i = 0. The sum of both dipoles dominates the structure 
of the distribution for transverse polarization of both quark and proton, as is seen in Fig. whereas 
the quadrupole-type term s t {2b l W — b 2 5 % ^)S^ H'^, is less prominent in our numerical example. We note 
that the two dipoles terms S l e lJ b 3 E' and s l e 13 V (E' T + 2H' T ) tend to cancel if quark and proton spin 
are opposite to each other, and the resulting density (not shown here) is rather sensitive to the precise 
values in the form factor parameterizations. Figure El finally shows the lowest moment of the vector 
field Fj,(x,b) describing the transverse quark polarization in a proton with or without transverse 
polarization. 

3 The spin matrix and positivity constraints 

In the previous section we have discussed the density of quarks with transverse or longitudinal polar- 
ization in a transversely or longitudinally polarized proton. Densities for arbitrary polarization states 



11 




Figure 4: Illustration of the lowest moment 4 j\ dx [F(x, b) + s l F^(x, b)] for n-quarks in a proton 
with transverse spin vector S = (1,0). The transverse quark spin vector is s = (1,0) in the left plot 
and s = (0, 1) in the right plot. 
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Figure 5: Illustration of the lowest moment j\ dx F^, (x, b) of the vector field describing the transverse 
polarization of u-quarks in an unpolarized proton (left) and in a proton with transverse spin in the 
x-direction (right). 
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can be obtained from the spin matrix in the light-cone helicity basis 



dz 
4tt 



ixp^ z 



p + ,0,A' q(z 2 )T x >xq(zi) P + ,0,A 



(32) 



with zf = Z2 = 0, z\ = Z2 = b, and z-jf = —z^ = \z~ . Here A' and A denote light-cone helicities 
of the proton states, and definite light-cone helicities A' and A of the quark are projected out by the 
Dirac matrices (see e.g. jH]) 



r ++ = 7 +(i + 75 ), 

r— =7 + (l-7s), 



r_ + = -ia +i (i + 7 5 ) 
r + _ = «T +1 (i-7 5 ). 



(33) 



The corresponding labeling of helicities in a handbag graph is shown in Fig. EK- The matrix Moy'A')(AA) 
reads 



/ 



H + H 



m 



-ie- i< t>—(E' T + 2H'n 
m V 

2 ( HT ~^2 AbHT 



-ie'^—E' 
m 

H-H 

2e -^.^-H" 

m z 

-ie-^—(E' T + 2H' T \ 
m V / 



^(E' T + 2H' T ) 2(H T 



m 

2e ^^H" 
m 2 



H-H 

ie i( t>—E' 
m 



4m 2 



^ b 



ie^—(E' T + 2H' T 
m 



H + H 



(34) 



if we order the proton-quark helicity combinations as (AA) = (++), ( — h), (H — ), ( )• Here GPDs 

are given in impact parameter space with the notation specified after (JHJ), and the azimuthal angle 4> of 
b is defined after (jllj) . The quark density for an arbitrary polarization state of proton and quark can 
be written as \ Sa'A'AA ( c A'A')*-^(A'A')(aa) c aa with coefficients normalized to Saa I c aa| 2 = 1- This 
implies that the matrix M^'A'XAA) must be positive semidefinite. Integration of (|34|) over b leads to the 
known spin matrix for the forward distribution functions fi(x), gi(x), h\{x) according to (|14j) . and the 
positivity of the corresponding eigenvalues gives immediately the Soffer bound 2|/ii(x)| < fi{x)+g\{x). 
Using the relations (|22|) we see that the matrix (|34l) of impact parameter dependent distributions is 
the exact analog of the spin matrix for transverse momentum dependent distribution functions, which 
was discussed in |2*5]. 

In order to simplify the following discussion, we change basis by multiplying (|34|) with the diagonal 



matrix D = diag(l, ie l 



-le 



-i<j> 



1) from the right and with from the left. This gives a matrix 



b_ 

m 



\ 



H + H 

m 

(yEj> + 2Hrp^j 
1 



m 
H - 

-2— g JEfy 



E' 

- H 



m 



(E' T + 2H' T ) 2 (H T --^A b H T 



2(HT-^A b H T ) ~{E' T + 2H 



T 



-2—H'Jp 
H-H 

m 



Am? 

I Erp + 2Hj 



-E' 



m 



H + H , 



(35) 



To compare with the matrices given in |28) one must take into account that in those papers the sign convention for 
the Sivers function f^ T is opposite to the convention from |21| used here, and that the rows of the matrices in those 
papers correspond to the indices (AA) rather than (A' A'). We thank A. Bacchetta for clarifying discussions on this issue. 
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(a) 



A 




A 



(b) 



Figure 6: Labeling of the proton and parton helicities in the matrix elements Mf/^x'VAX) m the quark 
region x > (a) and the antiquark region x < (b). 



which is purely real and depends on b but no longer on (p. Positivity of the upper left 2x2 sub-matrix 
of (|35j) leads to the bound 

ft,.,,, (36) 



— \E'\ < \/H 2 -H 2 , 
m 



which has been discussed in [22]. Using the eigenvalues of the full matrix (|35j) . we can tighten these 
bounds by including the tensor GPDs. With the abbreviations 



a± 
b± 
c± 



H ± (H T 



+ 



1 



m 
1 



2&T 2m2 



A b H T ) , 



H ± ( Ht T;H'n 



m 



m 



(e' ±E' t ± 2H' t 



the four eigenvalues read 



a+ + Jb% + c± , 



b 2 _ + c 2 



(37) 



(38) 



We see that they are related pairwise by changing the sign of all chiral-odd distributions. This is 
tantamount to multiplying the matrix (j35|) with diag(l, 1, —1, —1) from the left and from the right, 
which does of course not change its eigenvalues. Positivity of the eigenvalues (|35|) gives the bounds 



< a± , c 2 ± < a±-b± = {a± + b±)(a± - b±) , 

which in particular imply < a± + b± and < a± — b±. They explicitely read 



H T + 



1 



m 



1 



2m 2 



< H 



and 



b ^(E>±E> T ±2H> T ) 2 < (H±H T ±±H^^L 2 A b H T ) 2 -(H±H TT ±H! r 



(39) 



(40) 



(41) 



In the phenomenological study jHUl it was found that the bound (|3fij) can indeed be very restrictive 
(and thus helpful) in reconstructing generalized parton distributions from experimental data. The 
tighter bounds (|4()j) and (|41|) may therefore be of practical value even with limited information on 
the three additional chiral-odd distributions they contain. As in the case of the usual parton dis- 
tributions, renormalization of the operators in (|32|) may destroy the density interpretation of the 
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impact parameter distributions. Closer analysis reveals that the bounds following from positivity of 
the matrix AT(a'A')(AA) should be valid at a sufficiently high renormalization scale fi 2 [3^. They are 
stable under leading order evolution to higher scales, as shown in |32| . 

Since both experimental information and results from lattice QCD calculations are in the first 
instance given as a function of t, it is useful to have bounds also directly in momentum space. This 
can be achieved by applying to ((11)1 the method used in for the simpler bound (J3SJ), where the 
main ingredient is the Schwarz inequality. A method leading to the same results is to multiply 1)35)1 
from the left and the right with diag(l, mb, mb, 1). In the resulting matrix M only even powers of b 
appear. Integrating over b as 2tt Jq°° dbb M(b 2 ) = J d 2 b M(b 2 ), one obtains GPDs and their derivatives 
at zero momentum transfer t = 0. The result is still a positive semidefinite matrix, whose eigenvalues 
have the form (|38|) with 



b± 
c± 

This leads to the bounds 



,0 



H + H + Am —(H — H) ± 2(H T - 2H T ) 
at 



H + H — Am —(H — H)± 2(H T + 2H T ) 
at 



t=o 



2(E ±E T ± 2H T ) 



4=0 



( K E±E T ± 2H T 



< 



4=0 



(H + H± 2H t \ (Am 2 — (H - H) T 



(42) 



(43) 



4=0 



where both expressions in large parentheses on the right-hand side must be positive or zero according 
to our remark after ()39(). The condition < [H + H ± 2Ht\i=o is j us t t ne Soffer bound. 

Alternatively, we can multiply the matrix 1)35)1 from the left and the right with diag(m&, 1, l,mb) 
and then integrate over b as described above. The eigenvalues have again the structure of (|38|). and 
we obtain bounds 



(e±E t ±2H t Y < (H-H^-^r [ dt'H T (t')) (Am 2 ^-(H + H±2H T )^2H T 
\t=o l 2 00 



(44) 



4=0 



with both expressions in large parentheses on the right-hand side positive or zero. They contain a 
function which is non-local in momentum space, namely 



dtH T (t) = ATrH T (b = 0), 



(45) 



which can be traced back to the integral J d 2 b b 2 H^(b 2 ) in the derivation. 

For reasons which will become clear shortly, certain applications require bounds which do not 
involve the distribution H . One such bound is simply obtained by omitting the last term in 1)411) . 
which together with (I4()[l leads to 



m 



E' it Erp it 2H'rp 



< H± 



+ 



m 1 - 



T 



2m 2 



A b H T 



(46) 



To obtain a bound in momentum space we multiply ()46j) with mb, integrate over b, and use the 
Schwarz inequality in the forms / d 2 b f < J d 2 b\f \ and / d 2 bbg < (J d 2 bg) 1 / 2 (/ d 2 bb 2 g) 1 ^ 2 , following 
the method of This leads to 



(E±E T ±2H T y < (^H±H T ^^J t dt' H T {t')) (Am 2 ^(H±H T )T3H T 



(47) 



4=0 
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where both terms in large parentheses on the right-hand side must be positive or zero. Alternatively 
one can add to (|3lfl) the positive semidefinite matrix 



V 



H-H 


m 





—E 1 


H + H 





m 












H + H —E' 

m 








—E' H-H 

m 



(48) 



/ 



and proceed as above. One then obtains bounds analogous to (|41jl and to Q43|) and (|44[1 by the 
replacements H — > 2H, E — * 2E, H — > 0. They explicitly read 



(49) 



and 



[2E ±E T ± 2H T 
(2E + E T + 2H T 



J t=o 



t=0 



< 



< 



I ( /"/ ± Jl-r ) [^m z —H + 2H T 
1 /"* 



J t=o 



(50) 



t=o 



So far we have considered quark distributions. For antiquarks we define the matrix M^/x')(ax) ( x > b) 
with x < as in 1)32(1 . but with the Dirac matrices now reading 



r++ = -7 + (i-7 5 ), 
r_ = -7 + (i + 7 5 ), 



r_+ = i C 7+ 1 (i + 7 5 ), 
r+- = -m +1 (i-75) 



(51) 



instead of (|33|) . A global minus sign compared with the quark case arises because the order of the 
operators q and q in (|32[) has to be reversed to obtain a density operator for antiquarks. To understand 
the signs in front of 75, recall that antiquarks with positive helicity have negative chirality. One must 
finally keep in mind that the helicity index A refers to the parton on the left-hand side of the handbag 
diagram as shown in Fig. For antiquarks this parton is annihilated by the operator q, and not 
by q as for quarks. Comparing (jBTj) with (JHH]) . we find that the spin matrix M(a/a')(AA) (x, b) in the 
antiquark region x < reads as in (|34|). but with the signs of all GPDs except H reversed. It is 
positive semidefinite, and one readily obtains bounds for antiquarks analogous to those we have given 
for quarks. 

Let us finally address the question of positivity bounds for Mellin moments of generalized parton 
distributions at £ = 0, which are for instance relevant in lattice QCD calculations. As discussed 
in the previous section, these moments involve both the quark and antiquark regions, x > and 
x < 0. Clearly, it is only the sum of quark and antiquark densities and not their difference for 
which positivity is ensured. This leads us to consider the moments j\ dx x n_1 / (x , b 2 ) with even 
n for all distributions / except H, which is why we have derived bounds without H. To derive 
bounds for the even moments, we can add the positive semidefinite matrices Jq 1 dx x n ~ 1 M(x,b) and 
J -1 dx (— x) n_1 M(x, b). The result involves Mellin moments f_ 1 dxx n ~ 1 f(x,b 2 ) for all distributions 
/ except H, where instead one has dx x n_1 sgn(x)H(x, b 2 ), which is the matrix element of a 
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highly nonlocal operator. Positivity bounds are then obtained exactly as above, and we find that 
the inequalities (|46j). (|4"7|) and ((ST))) also hold when we replace all distributions with their even 
Mellin moments. 

It may be interesting to see whether the bounds which do involve H also hold when we replace this 
distribution by its Mellin moment j\ dx x n ~ 1 H(x, b 2 ), thus taking the "wrong" sign in the antiquark 
region x < 0, or whether the bounds given in this section also hold for odd Mellin moments. This 
would signal that the antiquark contribution to the moments in question is sufficiently small to not 
destroy positivity of the quark contribution, i.e. of the matrix Jq dx x n ~ 1 M(x, b). 

4 Equations of motion and distributions of twist three 

At the end of Sect. [2] we have encountered the total derivatives of the chiral-odd operators which 
define transversity distributions through the matrix element Fj,. Using the Dirac equation for the 
quark field operator, we can rewrite the local operators appearing in 1)25)1 as 

n 

= -2 q(iD + ) n l5 q + £ q{iD + r 2 a + m75 gG^iW+f^q , (52) 

i=2 

n 

= -2 q(iD + ) n q + £ q(iD+y- 2 a + ^ gG^iiD^q 

i=2 

+ 2m q qj + (iD + ) n - 1 q (53) 

for n > 1, where we have used d^[qa + ^ . . . q] = —iq(p-y + + 2D + — "f + P ) . . . q and [D^,D + ] = 
—igG^ + . Apart from the term proportional to the quark mass m q , the operators on the right- 
hand sides are of twist three. They are obtained by inserting covariant derivatives iD + into the 
pseudoscalar or scalar quark current and into the quark-antiquark-gluon operators qa + ^75 gG^ + q or 
qcr + n gG^ + q (which can be written in a number of ways using the relations <t A/x 75 = — ^e AMQ ^cr Q/ g and 
Q\fj, _ These quark-antiquark-gluon operators are chiral-odd partners of the operators 

obtained by inserting covariant derivatives into q^y + gG^ + q and <77 + 7s gG^ + q, which appear in the 
virtual Compton amplitude at twist-three accuracy and are well-known from inclusive deep inelastic 
scattering and from deeply virtual Compton scattering, see e.g. The forward matrix elements 

of the operators in (J53|) appear for instance in Drell-Yan pair production and have been studied in 
detail in |34| . A review of their properties, including their renormalization group evolution, can be 
found in Note that the derivative operator on the left-hand side of (|53|) does not contribute to 
forward matrix elements. 

The powers of covariant derivatives in (|52|) and ()53|) can be resummed to obtain nonlocal operators 
on the light-cone, which may be written as 

/oo r\?~ 1 

j*P+*- / dy - q { -y )W [-lz,y]T l ,gG^(y)W[y,^z] q ^z), (54) 

where z + = y + = and z = y = 0, the Wilson lines W are along light-like paths, and T and 
denote Dirac matrices. Local operators at position z = are then obtained by 

(P + ) n J dxx n - l 2 {x) = q{W + ) n - l Tq 

(P + ) n / dxx n - l Oz{x) = ^g(iB+) l - 2 r / , 5 G^ + (iS + ) n - i g (55) 

•* 8=2 



d u 



qa^iiD+r^q 



0„ 



qa 



+fl (iD 



+\n-l. 
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for n > 1. Note that the integral / dx 0%(x) is zero. The matrix elements of the nonlocal operators 
2 {x) and 03 (x) between nucleon states are parameterized by suitable generalized parton distribu- 
tions. Taking instead matrix elements between the vacuum and a meson state one obtains meson 
distribution amplitudes, and nonlocal versions of the equations of motion in (|52|) and (|53|) have been 
extensively used in this context jHSllHZJ- We note that the operators 03 (x) with = a +fl or <r +M 75 
involve only "good" components of the quark and gluon fields in the parlance of light-cone quantiza- 
tion and hence admit an interpretation in terms of parton degrees of freedom, unlike the operators 
2 (x) with r = 1 or 75, which are products of one "good" and one "bad" field component |341 138j . 

With possible applications to lattice QCD calculations in mind, we prefer here to work with the 
local operators in (|52j) and (|53j) and the form factors parameterizing the Mellin moments of GPDs. 
Instead of transforming (|25f) back from impact parameter to momentum space, we can directly use 
translation invariance to obtain (p'\d^O\p) = iA^(p'\0\p) for a local operator 0. From (|52j) and (|53|) 
we then obtain relations between the form factors of twist-two and twist-three operators. We give 
results for n = 1 and n = 2, their generalization to higher moments is straightforward. In contrast to 
the previous sections, we consider the general case where £ need not be zero. Using the constraints 
from parity and time reversal invariance, the quark tensor current can be parameterized by 

,Xfia/3 a p Xfia/3 a 

(p'\qcr X ^q\p} = ua x ^ 5 uA T10 (t)+u JLA u A TW (t)+u -^M u B T10 (t) , (56) 

Wr Zm 

where we use the notation of jSHj- An analog for the operator qa XfM q is readily obtained using 0^75 = 
— \ie^ a ^ a a 0. For the operator with one covariant derivative we have 

A^S, w (p'\qa X » Tj iD»q\p) = A A/1 S^ [^75 uP v A T20 (t) +u -^J- u P v A T20 (t) 

+ u eXfia(3 n Aalf3 u P»B T20 (t) + PalP u A u B T2l (t)] , (57) 

Zm m J 

where denotes antisymmetrization in A and /i and denotes symmetrization and subtraction 
of the trace. Comparison with readily gives 

J dxx n - 1 H T (x,£,t)=A Tn0 (t), J^dxx^HTix,^) = A Tn0 (t), 

J_ dxx n - l E T {x,i,t) = B Tn0 (t), J^dxxE T {x,i,t) = -2£B T21 (t) (58) 



for n = 1,2. In the forward limit we have Ht(x, 0,0) = hi(x), so that Arno(0) = f_ 1 dxx n 1 h\(x) is 
a moment of the usual transversity distribution. The contractions 

A^p'lqa+^D+r-'qlp) , A^qa+^D+^q^p) (59) 

needed for the equation of motion constraints (|5^)) and (|S1?)) respectively project out the form factors 

AtuU Btui and Axni, Ax n i, BTni- 

For the twist-two operators constructed from the quark vector current we have 

ia^A 

{p'\q^q\p) = u^uA w {t)+u— -uB 10 {t), 

Zm 

Sr,tf\<rfiD v q\p) = sJu^uP»A 20 (t) + u—^uP"B 20 (t)+u uC 2 (t)} (60) 
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and 



dx H(x, £, t) 



dxE(x,£,t) = Bio(t), 



^ dxxH(x,£,t) = A 20 (t) + 4^C 2 (t), J dxxE(x,Z,t) = B 20 (t)-ifC 2 (t). (61) 

Note that A±o(t) and Bio(t) simply are the contributions of the relevant quark flavor to the usual Dirac 
and Pauli form factors, respectively. In the forward limit ^4„o(0) = j\ dx x n ~ l f\{x) is a moment of 
the unpolarized parton distribution. 

We further parameterize the twist-three operators constructed from the quark scalar and pseu- 
doscalar currents as 

{p'\qiD^q\p) = mvrfuAsvi + \uia^ a A a u B S10 , 



S^{p'\qiD»iD»q[ 
(p'\qiD^ 5 q\ 
S^(p'\qiD^iD^ 5 q[ 



J fJ,U 



m vrfu P u A S20 + \uia^ a A a u P u B S20 + uu A^A V C S2 



h-un/suP" Bpio 



A u A P2l + \u lb u(p^P v B PW + A"A"Bp22 



(62) 



where we have omitted the argument t of the form factors for brevity. Following j!4l I39j we have 
assigned the subscripts of form factors such that the first subscript gives the spin of the operator (i.e. 
the number of Lorentz indices in the symmetrization and subtraction of traces) . The second subscript 
counts the number of factors A in the form factor decomposition whose Lorentz index corresponds to 
a covariant derivative on the operator side. In the forward limit p = p' only the form factors As n o(t) 
survive in (j62fl . They are the moments of the chiral-odd parton distribution e(x) defined in |34[ I2()j. 
given by As n o(0) = j\ dx x n e{x). Note that the local current qq without a covariant derivative 
(whose forward matrix element is related to the pion-nucleon sigma-term) does not appear in the 
constraints (|53j) . In other words, the equation of motion constraint involves xe(x) when resummed to 
x space, and thus is not affected by the S(x) singularity of e(x), discussed e.g. in the recent review [477] . 
If we finally define form factors for the quark-antiquark-gluon matrix elements as 



S^ip'lqa^gG^ql 
S^(p'\q^ a 75gG au q[ 



J flV 



J fJ,U 



2m fi/u P v A G2 q + uia^AaU P u B G20 + 2uu A»A U C G2 
2m vff<KU A u A G21 + n 75 n (pfp u B G20 + A^A^ B G22 



(63) 



the equations of motion embodied in (|52*|) and (|53|) give relations 



Asio 
Bsw 
Bpio 



A 



m 

m, 



10 



4m 2 

— - -Bio — f^-rio 
m \ 

2Atiq 



(Btiq + 2At\q ) , 
t 



2m? 



Atiq 



(64) 



for the lowest moments. At this level the quark-antiquark-gluon operators do not appear yet. 
In particular, at t = the first relation in (|64[) gives the well-known sum rule J_ 1 dx xe(x) = 
(m q /m) j\dx fi(x), where the integral over fi(x) at the right-hand side is just the number of valence 
quarks with appropriate flavor [HI]. The operators involving gluons do appear in the relations between 
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the second moments, 

Tfi ~b / ~ 

As20 — Aq20 = — - ^-20 — 1 — o ( Bt20 + 2At20 

Bs20 - B G 20 = — - B 2 o - (a T 2o - - — ~At2o) , Cs2 - Cq2 = — - C 2 , 

m \ 2m z / m 

Ap21 — Ag21 = ^2 #T21 , 

-Bp20 — Bg20 = 2At20 , -E>P22 ~ -E>G22 = ~Bt2\ ■ (65) 

In forward limit t = no connection is obtained between moments of the twist-three distribution 
e(x) and moments of the transversity distribution h\{x). Rather, form factors of the twist-three 
operators which survive in the forward limit are connected with form factors of the quark tensor 
current which decouple in that limit, and vice versa. Preliminary results on A^nO, AxnQ an d -Brno 
from lattice calculations |27| suggest that these form factors are rather large. If confirmed, this would 
imply that the twist-three combination As n o — AquO is quark mass suppressed at t = but no longer 
small for —t ~ m 2 , and that the form factor combinations BsnO ~ Bc n o an d Bp n Q — Bc n o are already 
large at t = 0. In other words, away from t = chiral-odd twist-three matrix elements would not 
generically be small compared with twist-two matrix elements. 

The relations (|fi4j) and 1)65(1 may be of practical use in lattice QCD calculations. Note that the 
form factors on the left-hand sides belong to operators with one covariant derivative more than 
those on the right-hand sides (counting the gluon field strength as the commutator of two covariant 
derivatives). Operators with more derivatives are less localized on the lattice and thus more affected 
by errors. The form factors on the right-hand sides of ([64(1 and 1)65(1 have been or are being calculated 
in lattice QCD. Together with lattice determinations of the renormalized quark masses (see e.g. |41j 
and references therein) one may thus use the equation of motion constraints to determine the twist- 
three form factors in (|64|) and the twist-three form factor combinations in l)65jl . Alternatively, one may 
evaluate the twist-three matrix elements on the lattice and use (|64j) and 1)651) as constraints to reduce 
the errors in the extracted form factors. Note that a separate determination of Ac-ni, Bc n i, Acni, 
Bcni would allow one to check the often-used Wandzura-Wilczek approximation, which assumes that 
matrix elements of (chiral-even or chiral-odd) quark-antiquark-gluon operators are small. 



5 Summary 

Generalized transversity distributions at zero skewness £ describe the density of transversely polarized 
quarks in the impact parameter plane. We have derived the corresponding expression (jBJ) and analyzed 
its detailed structure. The momentum space distributions Ht, Et and Ht at t = and £ = 
describe simple average features of this density according to (fT3|) . (fTTj) . (|T9|) . The formulae for 
the impact parameter density of quarks closely resemble those for transverse momentum dependent 
distributions. This resemblance exhibits for instance a correspondence between the Sivers function 
fyj, and the nucleon helicity-flip distribution E, and between the Boer-Mulders function hi and 
Et + 2Ht- It would be interesting to investigate the correspondence between impact parameter and 
transverse momentum dependent distributions at a dynamical level, as has been done for fyj, and 
E in |231 124j . The distribution of quarks in the impact parameter plane is no longer rotationally 
symmetric as soon as either the proton or the quark are transversely polarized, and the preferred 
direction of transverse quark polarization is not isotropic even in an unpolarized proton. Preliminary 
results of lattice QCD calculations suggest that such effects may be quite large. 

The impact parameter density of quarks for arbitrary polarization can be obtained from the spin 
matrix (|34j) . This matrix is positive semidefinite, which leads to simple bounds on generalized parton 
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distributions, as special cases of the general results in [32] • The most stringent inequalities hold in 
impact parameter space. A combination of Ht and Ht is bounded by H according to (|4Uj). and (|4Tj) 
extends the bound (j3(ij) previously given by Burkardt |29| . The size of the chiral-odd distributions thus 
has consequences also in the purely chiral-even sector, since it restricts the possibilities to saturate 
the inequality (|36j) . which involves only E, H and H. By suitable integration over the impact 
parameter, one obtains bounds in momentum space. Bounds can also be given for Mellin moments 
that correspond to the sum of quark and antiquark distributions. Since the axial current has different 
charge conjugation parity than the vector and tensor currents, this requires bounds without the quark 
helicity distribution H, like (jl6|) . (|4"Tj) and (fl9|) . (|50|) . Such bounds can for instance be applied to the 
results of lattice QCD calculations. It will be interesting to see by how much bounds are violated for 
Mellin moments corresponding to the difference of quark and antiquark distributions, since this is a 
measure for the importance of antiquark contributions in these moments. 

The divergence and the curl of the vector field F T (x,b), which describes transverse quark po- 
larization in the impact parameter plane, are matrix elements of the total derivatives of twist-two 
quark-antiquark operators. These derivative operators are related to twist-three operators via the 
QCD equations of motion, namely to scalar or pseudoscalar quark currents and to quark-antiquark- 
gluon operators. Such relations have been investigated for forward parton distributions and for meson 
distribution amplitudes in the literature. In (j64j) and (|65j) we give the corresponding relations for the 
form factors parameterizing the first two Mellin moments of generalized parton distributions. This 
can easily be extended to higher moments. Such relations may be of use for exploring the twist-three 
sector in lattice QCD calculations. 
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